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VO i ABSTRACT 
O 

^ ■ We construct models describing the velocity field in the infall regions of 

clusters of galaxies. In all models the velocity field is the superposition of a 

OS • 

radial systematic component which is assumed spherically symmetric, and a 
"noise" component of random nature. The latter accounts for the combined 
effects of small-scale substructure and observational errors. The effect of the 
+l5 . noise term is to smear out the caustic envelopes searched for by previous groups, 

resulting in models which resemble the observations quite well. When the 
. systematic component is known, the infall velocity and the mass profile of the 

^ . infall region of the clusters can be determined. 

Given a particular model, it is possible to calculate the distribution function, 
f(cz,9,m), at all points in the observable (cz, 9, m)-phase-space outside the 
virialized region. It is demonstrated, on simulated data, that the use of a 
maximum-likelihood estimator enables identification of the correct model, even 
at realistic noise levels. To minimize systematic effects due to deviations from 
spherical infall, observations from a number of clusters should be used when 
applying the method to a real dataset. Combining data from different clusters 
is a straight-forward procedure. 

In the models the mass-to-light ratio is allowed to vary freely with radius. 
Furthermore, the models do not require any knowledge of the virial region. 
Rather, such knowledge can be used to test the models ability to produce 
realistic results by comparing the mass estimated by the models at small radii 
to the mass estimates obtained for the inner region using the virial theorem. 
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As a corollary we show how distance information of the quality obtained 
using the Tully-Fisher relation enhances the likelihood signal. Such distance 
information might be used in the future to either strengthen the results for this 
type of model or allow more advanced models to be used (e.g. models breaking 
the assumption of spherical symmetry). 

The method relies on observations of redshifts of galaxies at large angular 
separations from the centers of the clusters. Currently such data exist only for 
Coma, to which we apply the method as an example. The estimated mass at 
r = 2.5h~ 1 Mpc falls in the range 1.1 — 2.4 x 1O 15 /i _1 M , depending on the model 
and the level of the noise term. We find r tU rn to be of the order lOh^Mpc, a 
result which is rather robust against changes in the model and subdivision of the 
data in various ways. We find 0.6 < < 0.8, where b is the bias parameter, 
when using the luminosity profile adopted in the models. Interpreting this result 
one must recall that the luminosity function of Coma is not well known in the 
infall region, and that ideally the method should be used on several clusters 
simultaneously to minimize the effects of deviations from the spherical model. 

The method provides a promising way of measuring the infall velocities in the 
outer regions of clusters. However, a robust determination of cluster properties 
must await future redshift observations of galaxies in the outer regions of a 
number of clusters. 

Subject headings: galaxies: clustering — clusters: mass, mass-to-light ratios 
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1. INTRODUCTION 

Current estimates of flo based on the mass-to-light ratio found for the virialized central 
regions of clusters of galaxies favour a low value (Carlberg et al. 1997a, David, Jones, & 
Forman 1995, Tyson & Fisher 1995, Squires et al. 1996, Smail et al. 1997, Fisher & Tyson 
1997), of the order f2 ~ 0.25. Estimates from much larger scales often favour higher values, 
including Q = 1, but the uncertainties are in general rather large (see Dekel, Burstein, & 
White 1996 for a recent review). 

If the global value of M/L is significantly higher than the value found for the central 
parts of clusters it would not be unnatural for the change to manifest itself in the infall 
region of the clusters. If, for example, the difference is due to environmental effects resulting 
in the baryons being "lit up" depending on the properties of the surroundings, one can 
argue that the infalling galaxies formed well away from the clusters and have yet to be 
strongly affected by cluster imposed tidal fields, hot cluster gas, etc. Furthermore numerical 
simulations indicate that the extent of the mass density profiles of clusters is related to f2 , 
more extensive profiles corresponding to higher values of Qq due to ongoing infall (Crone, 
Evrard, & Richstone 1994). The steepness is also related to the precise form of the initial 
spectrum of density fluctuations, however. 

Obtaining mass profiles of the outer regions of clusters is therefore highly relevant. 

Three main measures are being used for determination of the mass-to-light ratio of 
clusters: the velocity dispersion, the temperature of the hot X-ray emitting gas in the 
clusters, and the distortion of the shapes of background galaxies. The first two rely on the 
assumption of virial equilibrium and are consequently not the ideal tools for a determination 
of the mass at significantly larger radii than the virial radius, at least not without some 
modification of the approach governed by insight from N-body simulations. Observation 
of weak gravitational lensing using new 8-10-m telescopes and large array CCDs is a very 
promising technique which potentially can probe the cluster potentials far into the infall 
region. However, at present the mass estimates found from the weak lensing method 
generally apply to the inner regions. 

In this paper we shall develop a second method; we attempt to measure the systematic 
part of the velocity field of the galaxies in the infall region. 

As the peculiar velocities, relative to those expected for a pure Hubble flow, will be 
larger near a cluster than in the general field one might expect it to be a less difficult task 
to map the peculiar velocity field of the galaxies surrounding a cluster rather than the 
velocity field on scales on which linear theory applies. However, even though the infall 
signal should be rather strong , projection effects and our inability to measure distances of 
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galaxies precisely makes the determination difficult. 

Kaiser (1987) applied the spherical infall model (e.g. Lilje & Lahav 1991) to an 
(adopted) spherical cluster density profile and found that it manifest itself as a characteristic 
diamond shaped distribution in the distribution of galaxies in the (cz, #)-plane due to 
caustics, (here cz is the redshift and 9 is the angular separation between a given galaxy and 
the cluster center). 

Regos & Geller (1989, RG) extended the work of Kaiser and showed that the shapes 
and sizes of the caustics are a characteristic of the mass enclosed. RG assumed light to 
trace mass, in which case Qq becomes the only free parameter of the model, as the observed 
luminosity profile can be deconvolved and used to deduce the relative mass over-density 
as a function of radius. Fitting their models by eye to observations of the Coma-cluster 
(and a few other clusters) RG were not able to put new constraints on Qq, however. The 
data available to them were limited, and what they had appeared to be more noisy than 
predicted by their models. The density gradients across the caustics seemed to be weaker 
than predicted by their model. 

Van Haarlem et al. (1993, vH) adopted the basic assumptions of RG and refined their 
work on the Coma cluster in three ways: Firstly, they introduced an objective measure of 
the change in density between the region just inside and just outside the position of the 
caustics predicted by a given model, allowing a more rigorous fitting procedure. Using 
Monte Carlo simulated artificial clusters the new measure was shown to work well. Secondly 
they extended the number of observations available by performing a small redshift survey 
around the Coma cluster for ^[l ; 3°], adding 98 new redshifts to the Coma sample. Thirdly, 
using scans of Schmidt plates, they determined the luminosity profile of Coma as well as is 
currently possible (assuming spherical symmetry). 

Despite this the result was still inconclusive, which they noted was due to the caustics 
being smeared out in the observational sample, whereas the model caustics were very sharp. 
vH attributed the smearing observed in real clusters to a non-spherical mass distribution, 
and they demonstrated that an ellipsoidal mass distribution will produce a distinct smearing 
of the infall in the (cz, 6 l )-plane. Later work (e.g. van Haarlem & van de Weygaert 1993), 
show the non spherical distribution of the infalling material (infall along filaments), and the 
event-like nature of part of the infall, to be important factors, however. 

In this paper we show that keeping the spherical approximation for the infall, but 
adding noise to the models, noise which is known to exist in all clusters, and which 
represents the combined effects of the velocities induced by small scale density fluctuations 
and observational errors, can be responsible for much of the smearing of the caustics. We 
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are able to produce artificial clusters with our models which appear quite like real clusters 
in the (cz, 6>)-plane. 

The organisation of the paper is as follows: Section 2 describes the models, the 
maximum likelihood procedure used to fit models to data, and tests of the procedure 
on simulated data. In Section 3 the models are applied to the Coma cluster as an 
example. Section 4 contains a discussion of our findings and assumptions, and Section 5 the 
conclusions. 



2. MODELS AND TESTS OF MODELS 

We will construct models which allow us to determine the infall velocity profile and 
mass profile of clusters of galaxies beyond the virial region, based upon observations of the 
position on the sky, the line-of-sight velocity, and the magnitudes of individual galaxies in 
and around the clusters. 

We shall attempt to perform this task by first constructing a series of simple models 
of the velocity field of the region beyond the distance of shell crossing, with each model 
having only a few free parameters to be determined. Secondly, the distribution function, 
f(cz,9,m), which expresses the number density of galaxies expected at any given point in 
the observed region of phase-space is calculated for each model. Finally, the best-fitting 
parameters are found using maximum-likelihood techniques. 

With the model-parameters in hand, estimates of the basic properties of the outer 
regions of a cluster, such as the mass inside a given radius, the mass-to-light ratio (M/L), 
and the turn-around radius can be estimated, and eventual conclusions about Qo can be 
drawn. 

In all of the following we define the "infall region" as the region of a cluster which 
is outside the part seriously effected by processes associated with virialization and shell 
crossing and inside the volume outside which neighbouring clusters have a significant impact 
on the peculiar velocities of the galaxies relative to those induced by the cluster itself. We 
shall apply our models only to the infall region, as here defined. It contains galaxies which 
are falling towards the cluster for the first time as well as more distant galaxies which are 
still moving away from the cluster. 

We shall assume that the velocity field in the infall region can be adequately described 
by a systematic component which is the one we want to determine plus a 'noise' component 
of random nature. The latter allows for the effects of substructure, non-radial motions, and 
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observational errors. 

Following RG, we will let the systematic component of the velocity field be given 
by a spherically symmetric radial flow. Many clusters show evidence for deviations from 
spherical flow, but combining data from a number of clusters we expect the effects of 
non-sphericity to be minimized out, resulting in a spherical cluster while adding some extra 
noise to the infall signal. 



2.1. Infall Velocity Versus Density in the Spherical Approximation 

The spherical approximation is based on the constraint that one is dealing with a 
spherically symmetric density fluctuation and the corresponding purely radial velocity 
fluctuation in an otherwise perfectly homogeneous and isotropic background universe. In 
the limit of no shell crossings (i.e. there is no mass transport through the shell over the 
time span considered) the evolution of the radius of a shell and the current time are related 
through the following equations: 

v a a 

-=|^/(^) and tif o = _0(,7), (1) 

where: 

f(v) = I 1 " COS?? ' I > ° ^d g( V ) = { 71 ~ Si ^' I > ° , (2) 
Jyn \ cosh 77-I, b<0 yKIJ \ sinh 6<0' KJ 

where a, b, are constants related to the initial properties of the shell and to the type of 
cosmology (e.g. Peebles 1980). Throughout this paper we use lOOhMpc km^ 1 s" 1 for the 
value of the Hubble constant H . 

It follows that the ratio of today's velocity, v(r), of a given shell at distance r from the 
center of the perturbation to the velocity difference expected over the same distance in a 
situation with pure Hubble flow is: 

v i g{v)f'{v) 



H r t H ' 

whereas the ratio of today's mean density inside the same shell, p(r), to that of the 
background universe is given by: 

P_ = 2 g 2 (r]) 

Po (t Q H yn pin)' 

(e.g. Silk 1977). Remembering that: 

1 ^co S h(A -1), fi <l 



(3) 



(4) 



g, HQ — l, 
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and 

P « = ^G (6) 
it is possible to determine the mass profile in the infall region of a cluster if one can 
somehow measure its velocity profile. When combined with the luminosity profile of the 
cluster, M/L can be determined. 

Unfortunately it is currently not possible to measure the line-of-sight velocity as a 
function of distance directly from which v(r) could be obtained. Distance estimates for 
galaxies are far from being precise enough. One is forced to work primarily with angular 
position plus redshift data, and consequently v(r) must be obtained in a much more indirect 
way. 

In this work we do so by constructing a model which predicts the distribution of cz as 
a function of 6 on the basis of assumed functional forms of v(r), and subsequently we fit the 
distribution to the data obtaining a best-fitting v(r). 



2.2. Models of the Systematic Component 

It is non-trivial to write down the mass or velocity profile expected for the infall region. 
Most of our current knowledge about the density and velocity profiles of clusters comes 
from the inner regions. This is true both for observational and theoretical work (including 
N-body simulations). Simply extrapolating those profiles out to the infall region might 
lead to biased results. In the absence of knowledge of the actual profile, we have found 
it important to use a number of different model profiles in order to make sure the results 
obtained are not just artifacts of our ignorance. 

The adopted profiles should contain as few free parameters as possible, given the 
currently meager amount of observational data, the likely effects of the deviations between 
our spherical approach and reality, and the level of "noise terms" which smear out the clean 
signal of v(r). 

In terms of the ratio of the radial velocity relative to that expected for a pure Hubble 
flow at the same distance, H(r) = v(r) / (H r) = 1 — v pecu i iar (r) / (H r) , the first two of our 
adopted velocity profiles can be written as: 

H = l-(-^ (7) 

' turn 

H = l-/3exp( —log/3), (8) 

f'tum 

with each of them containing two free parameters which must be estimated. 
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Two additional profiles were given by specifying the average mass-density profile as a 
function of cluster-centric radius: 

P_ = 1 + ( PM _^(J_yi (9 ) 

PO PO r turn 

P _ P{ r turn) ^ 1" 
PO PO r turn 

For future reference we label the above four models VI, VII, DI, and DII, respectively 

As we have seen, for a given value of fl and H there is a one-to-one correspondence 
between a velocity profile and a mass profile. Working directly with a velocity profile 
emphasizes the fact that determination of the infall pattern can be considered a purely 
empirical task. Furthermore, it can be done independently of Qq. However, for certain 
combinations of the fitting parameters, the assumed profiles correspond to unphysical 
situations. The local density inferred from the velocity profile may become negative at 
certain radii within the region considered. It is straight forward to identify the subsets of 
parameters which correspond to unphysical situations, but we shall not pursue the point 
any further in this paper. Here we present the method rather than a determination of real 
cluster parameters. 

The profiles VI, VII, and DI all approach a pure Hubble flow at large distances, whereas 
the profile DII does not. Profile DII corresponds to a situation in which the surroundings of 
the clusters have a lower mass-density than the universal average. However, as we consider 
only the region near the cluster, one should see the inclusion of profile DII rather as an 
attempt to check what will happen if the density is allowed to vary significantly around the 
turn-around radius, rather than be forced to go smoothly towards the background density 
at the turn-around radius. While the VI profile corresponds to a faster infall at smaller 
radii, as expected for a "pure" infall situation, the VII profile reaches a minimum at a 
certain radius and then approaches zero as the radius decreases. The latter is in accordance 
with the behaviour found in N-body simulations, where (more or less by definition) the 
radial velocity is zero in the virial core of the cluster and is followed by a transition zone 
which can stretch outside the virial radius (see Figure 5 of Crone et al. 1994). We expect 
profile VII to be less sensitive to a misplacement of the inner cutoff in the data sample, 
should the cutoff be done at too small an angular separation. On the other hand, one can 
obviously not use the part of the VII profile inside the minimum to determine a density 
profile for the cluster, as it becomes difficult to judge how smoothly the mass estimate 
produced using the infall profile joins the mass estimate of the inner, virialized part of the 
cluster using profile VII. 
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2.3. The Distribution Function Corresponding to Purely Systematic Infall 

Using Figure 1 to define our coordinates, an observer located at a distance R from the 
cluster would find the line-of-sight velocity of a galaxy located at (9, s) to be (RG): 

v los = cz = H R{cos9 + (-| - cos9)H}, (11) 

Let the spatial number-density of galaxies be n ga i(r), the number of galaxies found in 
a small portion of the (cz, #)-plane then becomes: 



where: 

J = 



dcz 89 

ds ds 

dcz 89 

oe ae 



n(cz, 9) dczd9 = n ga i(r)s 2 sin 9\3\ 1 dczd9, (12) 
H R [f s 2 dH H r 

— {{ R- COs6) lWr + H- R } (13) 
is the Jacobian of the transformation between the two coordinate systems. This result was 
presented by RG (to within a typographical error). Zero points in J~ x {cz,9) make the 
distribution n(cz, 9) have a very characteristic "horn" shape. Although n naturally has to 
be integratable over a finite volume of phase-space, the high-density regions around these 
caustic surfaces do contain enough material for the surfaces to be easily recognisable even 
in a modest dataset (Figure 2a). Such a profile was noted by Kaiser (1987) and RG, and 
its use as a promising way of determining the infall velocities and enclosed masses as a 
function of radius has been investigated (RG, vH). 

However, a number of effects complicate the situation making a direct identification 
of the "horn" almost impossible in a set of real observational data. In most clusters only 
a small number of galaxies have been observed within the infall region. A dedicated 
observational program can change this picture for the better, but clusters contain a finite 
number of galaxies and identifying the "horn" may demand combining data for a number 
of clusters, even if every galaxy in the infall region has a measured redshift. Secondly, the 
velocity field is probably not perfectly radial, even in the region dominated by infall, as 
envisaged by the various degrees of substructure on small and intermediate scales seen in 
many clusters. That will tend to smear out the caustics (vH), as will observational errors in 
the redshifts. In the next section we shall include a term in the models which incorporates 
the effects of small-scale substructure and observational errors. 



2.4. Inclusion of the Noise Term 



Within the hierarchical picture of structure formation the dynamical time scales for 
the evolution of small-scale structures (galaxies and groups) are generally shorter than for 
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larger scales (clusters); a difference increasing if the power spectrum describing the initial 
density fluctuations on those scales becomes steeper. When the dynamical time scales are 
sufficiently different, the internal dynamics of pairs and groups of galaxies in the infall 
region can be considered independent of the slower formation of the cluster. If so, the 
effect of small-scale substructure on the velocity field of the infall region of a cluster will 
be isotropic and independent of the position of the galaxies with respect to the cluster. 
We shall assume this to be the situation, in which case a rough estimate of the effects of 
small-scale structure can be obtained by folding the distribution with a Gaussian smearing 
function, G(x,a) = l/(y/2na) exp(— 0.5x 2 /a 2 ), leading to: 

/+oo 
n(cz, 9')G(cz — cz, a cz )dcz' (14) 
-oo 

/+oo 
n gal (r[9,s'})G(cz[9,s'}-cz,a cz )s' 2 sm9ds', (15) 
-oo 

the latter being the form most easily handled in numerical calculations for relevant values 
of a cz . 

The observational data consist of inhomogeneous datasets obtained using different 
selection functions. In general, however, the selection functions are much nearer to being 
magnitude limited rather than volume limited, and we have to include the depth effect, at 
least for clusters as near as Coma. 

Let til(M) be the luminosity distribution of galaxies, and assume rii is independent 
of the position (within the infall region) relative to the cluster position. If the sample is 
magnitude limited, we can write the proper distribution function as: 

/+oo 
n ga i{r[9, s'])G(cz[9, s'] - cz, a cz )n L (M[m, s'])s' sm9ds' (16) 
-co 

and M(m, s) = m — 25 — 5log w s, with m and M being apparent and absolute magnitudes, 
respectively. 

This is our final distribution function. We shall now test its ability to infer correct 
parameters on simulated clusters, but first we demonstrate the effect of a cz . Figures 2a-2c 
show the distribution in the (cz, #)-plane with and without the noise term included. The 
improvement from adding the noise term is quite evident, when compared to a real case as 
in Figure 2b (Coma). Figure 3 shows the number density along a line-of-sight for various 
degrees of smearing. Figures 2a, 2c, and Figure 3 were constructed using profile VI, (3 = 1.4, 
and r turn = 9Mpc. 

Which value should be used for a cz ? There are not enough data to allow a cz to be a 
free parameter in the models. Rather, we find guidance from the observed pairwise velocity 
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dispersion (reduced to that of a single galaxy) at small scales. Unfortunately, there is no 
current consensus as to the proper level of the pairwise velocity dispersion. The recent 
determination by Markze et al. (1995) found &±-2 = 540 km s -1 , markedly higher than the 
previous value, o\-2 = 340 km s _1 , found by Davis & Peebles (1977). However, as shown by 
Markze et al. (1995), o"i_2 is strongly influenced by the fraction of highly clustered galaxies 
included in a sample, as these galaxies, despite being relatively few, have a significant effect 
through their large relative velocities. Removing all Abell clusters of richness class R > 1, 
they find (Xi_ 2 = 295 km s -1 . This could argue for the use of a low value. However, the 
current neglect of non-spherical effects and the future extra noise from superposition of a 
number of clusters calls for a value somewhat above the estimate from pairwise velocity 
dispersions. Currently we bracket the problem by using the two values, a cz = 300 km s _1 
and a cz = 400 km s _1 . 

In the above calculations we have adopted n ga i/n = 1 + £ C9 , with £ Cfl = (r/r c )~ 7 being 
the cluster-galaxy correlation function and using 7 = 2.5 (Lilje & Efstathiou 1988), n is 
the mean number-density of galaxies in the Universe. Whereas this is a natural choice in 
cases where little is known about the number density profile of the cluster in question, other 
choices are of course possible, e.g. using the observed (deconvolved) profile - though one 
should notice that it is very difficult to determine the luminosity profile of a single cluster 
in the region considered here. We have experimented using different parameterizations of 
n gaiij')- Fortunately, it turns that the method is rather insensitive to these changes. 



2.5. Fitting a Model to the Data 

We fit the models to the data by maximising the expression: 



N 

lnL(model) = y^lnp(c^, Qj, m^model), (17) 

i=i 

where: 

p(cz, 9, m model = — — -^- f / - 18 

J J J n acz (cz f , 9\ m')dcz'd0'dm' 

is the probability of observing a galaxy at (cz,9,m), given a particular model and its 
parameters, and the integral is over the region of phase space from which the data are 
drawn. N is the total number of galaxies in the sample. 

This procedure is superior to envelope fitting of the caustics, as all data points 
contribute to the fit, thereby making the most of the available information. Due to the 
inclusion of the noise term, p(czi, 9i, model) is a rather gently changing function in cz, 
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and in practice it turns out not to be necessary to discretize p(cz, 9, m) in 9 or cz before 
calculating L, even for rather small data samples. One should check, however, that no 
galaxies of highly unusual brightness for a sample of a given size and distance are included 
when fitting, as such a galaxy can have a large impact on the results. Such a galaxy should 
either be excluded or p(cz, 9, m) has to be discretized with respect to magnitude. 

An important, valuable feature of the model is the fact that it is not necessary to 
correct for foreground and background objects, as they are included in the model. This is 
not the case for members of neighbouring clusters however, and care should therefore be 
taken to select data from a region of phase-space surrounding each cluster which is small 
enough that nearby clusters will only have a minor effect on the velocities. On the other 
hand, in the absence of neighbour problems, the method produces better determinations of 
the model parameters when the cuts in cz are sufficiently widespread to include the caustics 
at all selected 9, and the cuts in 9 are wide enough to cover the whole region between the 
radius of no shell crossing and r turn . In each case one has to weigh these arguments and 
find a working compromise. 

2.6. Testing the Procedure on Simulated Clusters 

Previously we demonstrated that the inclusion of a noise term in the spherical infall 
model and a few further reasonable assumptions produce simulated clusters which 'to the 
eye' mimic real clusters quite well (Figure 2). 

We shall now test our maximum-likelihood estimator, by using it on simulated clusters 
in order to evaluate whether it is capable of identifying the correct parameters at all 
(recalling that the previous attempts to identify caustics for real clusters failed). If so, what 
is its sensitivity and how big are the errors arising from assuming wrong model parameters 
for a given type of cluster? 

We have performed numerous tests on simulated clusters for which all model parameters 
were known. In most cases the simulated clusters contained 143 galaxies (sometimes 286) 
chosen randomly according to the distribution function p(cz, 9, m\model) (Equation 18) and 
having 2 < 9 < 8, 4930 km s" 1 < cz < 8930 km s" 1 , a cluster-centric velocity of 6930 km 
s_1 > leg — 2.5, and a magnitude cutoff at 15.5. These parameters may appear odd, but they 
correspond to those of our prime sample of real data (Coma). 

The results are presented in Table la. In the table, column 1 indicates the systematic 
component with which the simulated clusters were made, column 2 gives the velocity 
dispersion used in the smearing function of the simulated clusters, columns 3 and 4 give 
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the systematic component and the velocity dispersion assumed when doing the likelihood 
analysis, column 5 is the number of simulated clusters on which the analysis was done, 
column 6 gives the average value of r turn and its dispersion from the maximum-likelihood 
analysis, and column 7 gives the average value of (3 (or 7 for the density profile models) 
and its dispersion. All simulated clusters made using the profile VI were given parameters: 
(rturn, P) = (9, 1.4), using VII: (r turn , (3) = (9,8), and using DI: (r turn ^) = (9.5,2.5). 

First we note that the adopted maximum-likelihood technique is indeed able to return 
the proper parameters, even for a cz as high as 400 km s _1 . As might be expected the 
precision depends on a cz , the match between the model assumed for the simulated cluster 
(reality), and the model used for analysis. However, even in cases of mis- match or high a cz 
the results are still reasonable. 

Figure 4 gives an example of a contour plot for the maximum-likelihood levels found for 
a random simulated cluster (a cz = 300 km s _1 ). The plot shows ln(L), the cross identifies 
the position of maximum-likelihood, and the levels are separated by ln(2) and constantly 
decreasing as one gets farther from the peak. There is some degree of degeneracy between 
r tU rn and 7 (and similarly for velocity profile models between r turn and (3). This reflects 
the fact that for certain correlated changes of (r tU rn,l) the corresponding changes in the 
velocity profile (and density profile) are rather small. The uncertainty of an individual 
"measurement" in some cases is rather large, and ideally a number of clusters should be 
used when determining model parameters for a set of real data. 



2.7. Extension of the Model to Include Distance Information 

Obviously one cannot determine the parameters of an infall profile, v(r) without 
knowledge of both v and r. In the model above that information comes about rather 
indirectly in the form of distribution functions calculated on the basis on an assumed 
profile v(r) (and other model parameters) and the observed values of cz and 9. Due to 
the existence of multiple solutions for the distance for certain combinations of (cz, 9) even 
without noise, and due to the assumption of a random scatter in cz of the order a cz > 300 
km s -1 , the probability distribution, p(s\cz, 9, to), of galaxy being at a particular distance 
given the observed (cz, 9, to) is very broad. For most galaxies it is even impossible to tell 
whether they are in front of or behind the cluster. But as has been demonstrated it is still 
possible to reliably estimate the parameters of the infall profile. 

It seems worthwhile to test whether use of observed distances leads to an improved 
determination of the model parameters, even though the uncertainty in distance estimate 
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for galaxies is at present quite considerable relative to the sizes of the clusters in question. 
(If that was not the case the whole procedure presented in the previous sections clearly 
becomes obsolete, and one could do much better using another approach.) 

Therefore, assume that for a given distance estimator the probability that a galaxy 
of true magnitude M t and inferred magnitude M Q is given by G(M a — M t ,a m ). For the 
Tully-Fisher relation a rn is of the order 0.38 (c.f. Willick et al. 1997). With s being the true 
distance, the probability of observing a galaxy at distance s Q becomes: 

Q -y 25iog^ (s o /s) 

p(s \s) = = e . (19) 

Including this term in the expression for n(cz, 6, m) (Equation 16 and further) provides 
a simple way of including the distance information in the likelihood estimates. (As the 
above expression is already normalized and does not depend on model parameters other 
than a m , the extra integration over s Q in the normalization of n(cz,9,m, s ) is trivial as 
long as no cuts are made in s Q space when selecting the sample.) 

Table lb shows the result of including distance information in the tests on simulated 
clusters. Except for the added distance information (assuming a m = 0.38) for the simulated 
clusters and the inclusion of the distance term in the distribution function used in the 
maximum-likelihood analysis, everything is performed exactly as outlined in the previous 
sections, including the format of the table. Clearly distance information has a significant 
impact even for realistic values of a m . It would be of great value to have distances for the 
galaxies which are to be used in deriving infall profiles, despite the ~20% uncertainty of a 
single distance measurement. It is not clear whether the best observing strategy is to start 
measuring more distances to galaxies in the infall regions of a few clusters or to measure 
many more redshifts of galaxies in a larger number of clusters. 



2.8. Combining Data from Different Clusters or Datasets 

As long as the parameters to be determined can be assumed the same within different 
clusters or datasets the information can be combined by simply adding the likelihoods 
found for the individual datasets before searching the best fitting parameters. It is thus 
possible to combine, for example, a dataset with distance information with a larger one 
without distance information for a given cluster, even if they were made using different 
selection criteria or to combine different clusters observed using different selection criteria. 
Clearly the galaxies within a given dataset must have been chosen randomly from within 
the redshift, angular separation, and magnitude bins describing that dataset. 
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If the clusters in question are not "equal" , which in general will be the case, it becomes 
necessary to rescale them before performing the maximum-likelihood analysis. This can be 
done easily and with good results, as demonstrated by Carlberg et al. (1997b), who studied 
the virialized regions. 

3. RESULTS FOR COMA 

We shall now apply the models presented above to the Coma cluster. As mentioned 
above one should be cautious when using data for one cluster only. Unfortunately, Coma is 
currently the only cluster for which we have access to the type of data needed for a large 
part of the infall region. 

Since the model includes contributions from background galaxies, but not from other 
clusters, it is important to select a region of phase space which is not 'contaminated' by 
other clusters, yet contains as many galaxies belonging to the infall region of Coma and 
its outskirts as possible and simultaneously is large enough to constrain the sought-after 
parameters efficiently 

Figure 5 shows Coma viewed in the (cz, #)-plane. Filled circles are data from zcat 
(Huchra et al 1993) and open triangles are from vH. The rectangle inside the dashed lines 
defines the region of phase space we decided to model. It is a bit narrower in cz than we 
would like, but this was done in order to avoid including a significant number of galaxies 
from the 'filament' extending from the cluster A1367. 

The inner cutoff has to be large enough so that virialization and shell crossing have 
no effect on the galaxies selected. We use 9i ow = 2° which corresponds to r pa 2Ah~ 1 Mpc. 
The outer cutoff is 8°, large enough to contain most of Coma, yet small enough to avoid 
inclusion of galaxies heavily interacting with the surrounding clusters. The precise values 
of the two cutoffs are not important, and we have checked this by also running our models 
with cutoffs at 3° and 10°. 

It is well known that Coma is not perfectly spherical and that its early- and late-type 
galaxies are not distributed similarly. This has been demonstrated convincingly, for 
example, by Doi et al., (1995). To study the effects of this substructure in Coma on our 
analysis, we subdivided the dataset in three different ways: firstly into early- and late-type 
galaxies (according to Hubble type), secondly into inner and outer galaxies (according to 
6), and thirdly in high- and low-density regions, according the whether or not the galaxies 
are near or far from the two filaments (in 0) visible in Figure 6. In the latter two cases 
the cuts were made such that each splitting resulted in approximately the same number of 
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galaxies in the two parts. 

Despite going to much fainter luminosities, there is a distinct drop in the number of 
galaxies per magnitude bin at 15.5 in the full sample. Therefore, we have made a lower 
cutoff in luminosity at m=15.5. Furthermore, we decided to ignore the extra data included 
in the vH sample, as only few of them are found within the region of phase space considered 
here. 

The results are presented in Tables 2 and 3, while Figures 7a - 7c are examples of the 
likelihood contours found (with symbols as in Figure 4). In Table 2 the left half shows 
results for the DI model, while the right half are results for the DII model. The fits are for 
all 143 galaxies in the dataset. The value of f2 shown on the left is the one assumed in 
the analysis in order to convert a mass profile to a velocity profile. In Table 3 the left half 
contains results for the VI velocity profile, the right half for the VII velocity profile. The 
uppermost line gives the result for all 143 galaxies, whereas the results below are for the 
various subdivisions described above. 

We note how robust the estimates of the turn-around radius are to changes in the 
assumptions concerning Qq, o~ cz , and the velocity profile. For the total dataset the full range 
of values returned is 9.0/i _1 Mpc < r turn < 10.0/i _1 Mpc. Table 2 illustrates the effect of 
changing a cz . An increase in a cz will increase the width of the caustics for a given velocity 
profile, hence a fixed dataset results in a lower estimate of the systematic infall velocity at 
a given distance from the cluster. 

The only subsets of Coma data which result in a clear change in the velocity profile are 
the "high-density" galaxies and the "inner" galaxies. These are exactly the subsets that we 
would expect to deviate a priori from the analysis, the high density galaxies because local 
effects are likely stronger in the dense surroundings and because the marked deviation from 
our assumption of spherical symmetry (e.g. n ga i(r) being spherically symmetric). The inner 
galaxies deviate because the strongest signal in that region comes from the galaxies falling 
fastest towards the cluster, thus providing little information towards a determination of the 
turn-around radius. 



3.1. Implications of the Results for Coma 

What can we deduce from the Coma results? Firstly, we do not have enough data to 
arrive at any firm conclusions. The parameters found are not strongly constrained, and 
working on a single cluster our assumption of spheridicity might well break down. 
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Having said that it should be noted that the velocity profiles found are very similar 
for a given value of a cz . Not only is this true among models and cosmologies, interestingly 
enough it is also the case for the subdivided sets of Coma, with the exception of precisely 
those sets for which we would expect deviations. It seems our results are not completely 
arbitrary! 

As a further check we can compare the mass estimates provided by our models at 
a small radius to the virial estimate of the mass of the Coma cluster. This is done in 
Figure 8, where we plot the estimated masses at r = 2.5h~ 1 Mpc and r turn using the 
DI profile results. Circles are for a cz = 300 km s _1 , triangles for a cz = 400 km s _1 , 
and filled symbols for the smaller of the two radii. The inner mass estimates are nearly 
insensitive to the adopted value of f2 , as one would expect. The mass estimates at 
2.5h~ 1 Mpc are around 2 x 1O 15 /i _1 M for the low-velocity dispersion models and around 
1.5 x 10 15 /i~ 1 M o for the high-velocity dispersion models. Based on X-ray data Hughes 
(1989) found M(r < 2.5h~ 1 Mpc) = 0.5 - 1.5 x 1O 15 /i _1 M , with preferred values in the 
range M(r < 2.5h~ 1 Mpc) = 0.8 - 1.1 x 1O 15 /i _1 M , while Watt et al. (1992) estimate 
the total mass to be 1.3 x 10 15 /i _1 M Q . Using galaxy line-of-sight velocities, The & White 
(1986) find M(r < 2.7hT 1 Mpc) = 0.3 - 2.5 x 1O 15 ^ _1 M , with a value of 0.95 x 1O 15 /i _1 M 
assuming mass and galaxies are similarly distributed. Colless & Dunn (1996) find 
M viria i = 0.9 x 1O 15 /i _1 M . As emphasized by these and other authors, determinations 
based solely on galaxy velocities yield large uncertainties. Recall that our results do not rely 
in any way on data from the virialized region (except for providing the center of Coma). The 
good agreement between our results and the results quoted above is strong evidence that 
the method works, notwithstanding our "testing" laboratory with its simulated clusters. 
Performing a similar analysis, but using results found for slightly different cuts in angular 
selection (3-8°, 2-10°, 3-10°), we find M(r < 2.5h' 1 Mpc) = 1.1 - 2.4 x 1O 15 /i _1 M , again 
indicating that our results are far from accidental. 

If we had good knowledge of the luminosity profile of Coma, its M/L ratio as a function 
of radius could be derived for the whole infall region when combined with our results. 
Unfortunately, such a profile does not yet exist. vH did a very careful study of the light 
profile of Coma, using Ca and APM data separately to derive two light profiles, but they 
are based solely on data from the region 6 < 2°, making extrapolation to r tU m ~ 10/i _1 Mpc 
too uncertain. 

Alternatively, we attempt to use the luminosity profile, n ga i(r), specified in the 
models. Kent & Gunn (1982) estimated the total luminosity of Coma inside 9 < 3° to be 
L B = 4 x l0 12 h~ 2 L Q . This value can be used for normalization of our luminosity profile, 
and we can derive estimates of the total luminosity inside a given radius. As an example, 
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Figure 9 shows the mass-to-light ratio calculated this way at r = 2.5h~ 1 Mpc and at r turn 
(filled versus open symbols) for the DI profile for both the low and high value of the velocity 
dispersion parameter (circles and triangles). 

In the slightly non-linear regime one has: 

H r~ 3 (l + A)0-25' ^ 

where A is the relative mass over-density (Yahil 1985). That the luminosity does not 
necessarily follow the mass can be taken care of through a bias parameter: A^ = 6 A. In 
this case the peculiar velocity scales approximately as fl^ 6 /b ' 75 when A^ is substituted 
directly for A in Equation 20, in the outer infall region where the average over-density is 
larger than one, yet small enough that Equation 20 still provides a good fit to the real 
non-linear infall solution (see RG Figure 5). 

In Figure 10 we show some of the velocity profiles obtained using our model as well a 
set of curves based upon Equation 20 and using n ga i(r) to provide an estimate for A^. The 
thick solid lines correspond to predicted velocity profiles based on the light profile assuming 
tt° 6 /b - 75 = 1.0, 0.8, 0.6, 0.4, and 0.2, where 1.0 is the lowest line. The thin solid lines are 
from models VI and VII, the thin dashed lines the DI profile for both values of a cz , and the 
thin long-dashed lines are for the DI profile as above, but for 9 in the 3° - 10° interval. From 

Q0.6 

Figure 10 we see that 0.6 < ^5 < 0.8, but we caution the reader that the basis for the 
luminosity profile we use, the cluster-galaxy correlation function, is not necessarily a very 
precise description of the light profile of Coma throughout the infall region. While it serves 
in our models as a reasonable luminosity profile (since the models are not very sensitive 
to this), one should be careful in interpreting variations of the M/L ratio of Coma based 
on this profile. Also recall that our mass estimates for the inner region were somewhat too 
high, indicating the estimated infall velocities could be systematically a bit too large. 

Currently much work is being done on the relation between cluster profiles and 
cosmology (e.g. Cole & Lacey 1996, Navarro, Frenk, & White 1997). Nearly all the work, 
however, applies mainly to the virialized part of the clusters. It would be very useful to 
extend the work to the infall region, in which case the density profiles derived using models 
such as those presented here could be compared with results from N-body simulations. 

The values of r turn for Coma inferred from our models show little variance. Given 
some well-defined inner radius (e.g. the virial radius, r vir ) or the radius at a given average 
over-density, one might intuitively expect that the ratio r turn /r vir would allow discrimination 
between different cosmological models. There is an indication from Figure 5 of Crone et 
al. (1994) that this may be the case, but it is probably only through future high-resolution 
N-body simulations that one can hope to answer such a question. 
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4. DISCUSSION 

The main difference between the models presented here and previous attempts to 
determine the masses of clusters on the basis of their infall velocity profiles is the inclusion 
of a noise term, which represents the effects of small-scale structure and observational 
errors, and the use of the full distribution function, f(cz,9,m), rather than the caustic 
envelope, when fitting models to data. Furthermore, our modelling does not rely strongly 
on the light profile, which is not well determined in the infall region of individual clusters. 

Visual inspection reveals that in the limit of a no-noise term the model distribution 
and the real distributions are dissimilar (Figure 2). With hindsight it is therefore not 
surprising that RG and vH were unable to decide upon a best-fit model. Had they found 
one it might have resulted in false conclusions, as the effect of small-scale substructure and 
observational errors, represented through a cz in our models, is to increase the apparent 
width of the caustic surfaces. Without accounting for this effect the inferred systematic 
velocities become too high. 

With the noise term included in our models, the (cz, 9, m) distributions become 
acceptable to the eye (Figure 2), but as discussed earlier such a smearing of the caustics is 
also expected from an aspherical infall. Which of the two interpretations is correct? Most 
likely a combination, as both explanations are based upon known physical effects. It is well 
known that clusters are not perfectly spherical (see Figure 6 of the Coma cluster) which 
must result in a systematic flow which is not perfectly radial if the distribution of the dark 
matter is anywhere near the distribution of the galaxies. On the other hand, galaxies do 
induce substantial non-radial motions upon one another, even on scales so small that the 
effect clearly is not accounted for by a non-spherical cluster potential (e.g. in pairs and 
groups of galaxies) which should be combined with an observational error of the order ~100 
km s _1 . Furthermore, the level of o cz which leads to models resembling Coma is of the 
order 300-400 km s _1 , which is very close to the observed pairwise velocity dispersion for a 
single galaxy, which we expect to be a good measure of the "noise" on small scales. 

Are the effects of overall non-sphericity, be it a non-spherical potential resulting in 
a non-radial flow or a radial flow but mainly along a few filaments, so strong that they 
hamper the use of the spherical infall model, or the use of the particular models presented 
here? There is no doubt that using the spherical approximation for individual clusters one 
risks drawing false conclusions, as demonstrated by van Haarlem & van de Weygaert 1993 
on N-body data. They found everything from good to poor agreement between a spherical 
infall solution and the actual radial velocity distribution of the various clusters in their 
simulations. Thus to model individual clusters one will in general need to go beyond the 
spherical approximation. An attempt to do that has been made by Diaferio & Geller (1997). 
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But evidently an ensemble cluster, made by the superposition of a number of clusters, 
must be spherical. This is found to be the case for clusters in N-body models (van Kampen, 
1995). As one does not have spatial 3D information about real clusters, one might worry 
that selection effects makes an "observed" ensemble cluster non-spherical, but rather 
elongated along the line of sight. We believe that such problems can be minimized by 
using an X-ray selected cluster sample. Work is in progress, using the simulations by van 
Kampen (1995), to quantify to which degree there might be a problem if using an ensemble 
of clusters selected on the basis of "optical" data. 

The ensemble cluster being spherical, the relevant question is whether the systematic 
flow of the galaxies in the ensemble cluster can be described by the spherical approximation 
and related to cosmology. Lilje & Lahav (1991) found this to be the case. They compared 
the predictions of the spherical infall model applied to density peaks in a random Gaussian 
field to the outcome of N-body simulations and found the analytic models to give a good 
representation of the density and velocity profiles found in the simulations. 

We conclude that models like the ones presented here will work correctly when applied 
to a few clusters simultaneously. Further evidence comes from the mass estimates of the 
virial region of Coma produced with our infall method which are fairly close to the estimates 
produced based on observations of the inner region of Coma. 

How superposition of cluster data will effect the best value of a cz to use in the models 
is hard to say. On the one hand, one might expect that if non-radial flows are significant 
in individual clusters it will lead to an increase in the best value to use when considering 
a superposition. On the other hand, the values of a cz might be high already. Assume, for 
example, that deviations from spherical symmetry smear out the infall profile of Coma. 
Then the fact that our models represents Coma quite well could indicate that the level of 
a cz adopted is higher than what corresponds solely to small-scale substructure. Hopefully, 
future observational work on the level of the pairwise velocity dispersion on small scales 
can help answer such questions. It is an important parameter in the models, and it should 
ideally be determined without using the models themselves. For the purpose of these 
models one needs estimates of cr pairwise which excludes the central regions of clusters, as we 
need the velocity dispersions on small scales in the probably much "colder" infall regions. 

The best way to improve the models of the infall region is, we believe, by means 
of studying N-body simulations with emphasis on this region. As the typical systematic 
velocities in the infall region can easily be as small at the internal velocity dispersions 
of galactic dark matter halos particular care has to be taken the problem of identifying 
"galaxies" and measuring the velocity field of only those becomes critical. Simulations 
including simplified galaxy formation, like the ones done by van Kampen (1995), seem 
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an ideal basis for such a study. Potentially N-body models can provide knowledge about 
eventual systematic effects such as asphericity and about the way in which a cz might 
change due to superposition of clusters, or with the position relative to the cluster center. 
Also, given that we have found that the value of the turn-around radius of Coma is very 
insensitive to our modelling, it would be of great interest to identify in N-body simulations 
any systematic relation between the outer properties of clusters (such as the turn-around 
radius) and the inner properties (parameters which can be determined on the basis of virial 
estimates or similar) which relate to cosmology. 

5. CONCLUSIONS 

We have modelled the velocity field of the infall region of clusters of galaxies by a 
superposition of a systematic component, which is spherically symmetric and based on the 
spherical infall approximation, and a noise component, which is isotropic and independent 
of position in the cluster, which represent the combined effects of small-scale substructure 
and observational errors. These random velocities are assumed Gaussian and their level 
described by a cz . In the models the mass-to-light ratio is allowed to vary freely. Knowing 
the systematic component one can deduce the velocity and mass profile of the clusters. 

For a cz = 300 — 400 km s' 1 our models resemble Coma quite well. This level 
corresponds to the observed pair-wise velocity dispersion of galaxies on small scales. The 
effect of the noise term is to smear out the caustic surfaces searched for by previous groups. 
The models allow calculation of the phase-space density of galaxies expected at any point 
in the observable (cz, 8, m)-space outside the region of shell crossing. The likelihood of 
obtaining a given set of data from a particular model and its parameters can be determined. 
It is demonstrated on simulated data the maximum-likelihood procedure allows one to 
identify the correct model parameters, even for realistic noise levels and sample sizes. 

To avoid the effects of deviations from sphericity and to enhance the likelihood 
signal, data for a few clusters should be superimposed when applying the model to real 
observations. This is technically a simple process, but currently sufficient observational 
data are not available. 

The models do not in any way rely on data for the virial region of the clusters, other 
than providing the center of the clusters. Hence, the model determination of the cluster 
mass at small radii can be compared to the virial estimates as a test of the method. 

Finally, it is demonstrated that the inclusion of distance information of the precision 
which can be obtained using the Tully-Fisher relation strengthens the likelihood signal 
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significantly. Potentially this could be used to allow deviation from the assumption of 
spherical symmetry. 

We have applied our models to the only cluster for which we have enough data far out 
in the infall region, namely Coma. The estimated mass at r = 2.5h~ 1 Mpc falls in the range 
1.1 — 2.4 x 1O 15 /i~ 1 M , depending on the model and the level of the noise term. Higher a cz 
produces a lower mass estimate. We find the turn-around radius, r turn , to be of the order 
10h~ l Mpc, which is fairly insensitive to changes in the models, the noise level, and even to 
subdivision of the data. 

Using the cluster-galaxy correlation function to provide an estimate of the light profile, 

Q0.6 

we find for Coma 0.6 < < 0.8, b being the bias parameter: bA mass = A Ught , where A 
is the relative over-density. We hesitate to draw any firm conclusion on the basis of this 
result, as long as the present analysis was performed on data from only a single cluster. 

However, we have demonstrated that the method does indeed work, and it provides 
a promising way of measuring the infall velocities in the outer regions of clusters. More 
observations of redshifts in the outskirts of clusters are needed as well as a dedicated study 
of the properties of the infall region of clusters from N-body simulations. 
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Table la. Fits to Simulated Clusters. 



Simulated 
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n —1i 
[km s J 




n —It 
[km s J 






[MpcJ 




DI 


300 


DI 


300 


300 


143 


9.70±1.69 


2.53±0.52 


DI 


300 


DI 


300 


150 


286 


9.64±1.17 


2.50±0.37 


DI 


400 


DI 


400 


300 


143 


9.59±2.05 


2.41±0.61 


DI 


400 


DI 


400 


150 


286 


9.69±1.41 


2.45±0.41 


DI 


400 


DI 


300 


300 


143 


9.38±2.73 


2.67±0.73 


DI 


300 


DI 


400 


300 


143 


9.93±1.50 


2.27±0.47 


VI 


300 


VI 


300 


70 


143 


9.01±1.65 


1.44±0.33 


VI 


400 


VI 


400 


100 


143 


8.93±2.29 


1.39±0.35 


VII 


300 


VII 


300 


20 


143 


9.20±1.23 


8.31±2.51 


VI 


300 


VII 


300 


20 


143 


9.49±1.80 


7.51±3.00 


VI 


400 


VI 


300 


100 


143 


8.47±2.73 


1.57±0.38 


VI 


300 


VI 


400 


100 


143 


9.70±1.73 


1.10±0.37 



Table lb. Fits to Simulated Clusters, Including Distance Information. 
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[km s _1 ] 




[km s _1 ] 






[Mpc] 




DI 


300 


DI 


300 


300 


143 


9.51±0.91 


2.51±0.33 


DI 


300 


DI 


300 


150 


286 


9.52±0.60 


2.51±0.23 


DI 


400 


DI 


400 


300 


143 


9.33±1.04 


2.49±0.38 


DI 


400 


DI 


400 


150 


286 


9.41±0.69 


2.49±0.25 


DI 


400 


DI 


300 


300 


143 


9.43±1.26 


2.68±0.41 


DI 


300 


DI 


400 


300 


143 


9.37±0.83 


2.35±0.31 


VI 


300 


VI 


300 


70 


143 


8.99±0.81 


1.42±0.17 


VI 


400 


VI 


400 


100 


143 


8.78±1.17 


1.39±0.23 


VII 


300 


VII 


300 


20 


143 


8.89±0.99 


8.90±2.38 


VI 


300 


VII 


300 


20 


143 


10.40±0.87 


5.28±1.31 


VI 


400 


VI 


300 


100 


143 


8.95±1.44 


1.53±0.25 


VI 


300 


VI 


400 


100 


143 


8.56±0.86 


1.25±0.21 
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Table 2. Fits of Dl and D2 to Coma. 





0"cz = 


300 km s~ l 


&cz = 


400 km s" 1 


0~cz = 


300 km s~ l 


&CZ = 


400 km s" 1 


^0 


7 


^turn 


7 


Tturn 


7 




7 


f'turn 






[Mpc h- 1 ] 




[Mpc h- 1 ] 




[Mpc h- 1 ] 




[Mpc h- 1 ] 


0.01 


2.8 


9.5 


2.5 


10.0 


2.8 


9.5 


2.5 


10 


0.2 


2.6 


9.5 


2.4 


10.0 


2.8 


9.0 


2.5 


9.5 


0.4 


2.5 


9.5 


2.3 


10.0 


2.6 


9.0 


2.4 


9.5 


0.7 


2.5 


9.5 


2.2 


10.0 


2.7 


8.5 


2.4 


9.0 


1.0 


2.4 


9.5 


2.2 


10.0 


2.6 


8.5 


2.3 


9.0 



Table 3. Fits of VI and V2 to Coma. 



Component 


e 


N 


f'turn 


P 




P 




degrees 




[Mpc h- 1 ] 




[Mpc h- 1 ] 


[Mpc h- 1 } 


All 


2-8 


143 


9.5 


1.6 


10.0 


9.0 


Early types 


2-8 


51 


9.0 


1.7 


10.0 


8.0 


Late types 


2-8 


92 


10.0 


1.7 


10.5 


10.5 


High density 


2-8 


74 


13.5 


1.4 


14.0 


7.0 


Low density 


2-8 


69 


9.5 


1.4 


10.0 


6.5 


Inner region 


2-5 


73 


8.0 


2.0 


8.5 


12.0 


Outer region 


5-8 


70 


10.5 


2.0 


10.5 


12.0 
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Fig. 1. — A diagram showing how the various geometrical quantities in the models are 
defined. 
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Fig. 2. — (Q-.cz) diagrams of, (a) model based on profile VI with r turn — 9, (5 — 1.4 and no 
smearing term (Equation 12 but including magnitude term as in Equation 16), (b) Coma in 
similar diagram (only galaxies having apparent magnitude < 15.5), and (c) as (a) but now 
including a smearing term in cz of strength a cz = 300 km s _1 (Equation 16). 
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Fig. 3. — n(cz\9, a cz ) curves for 6 = 4° and different values of a c 
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Fig. 4. — Likelihood contours for a simulated cluster constructed using model DI with 
7 = 2.5, r turn = 9.5, and Qq = 0.4. The cluster contains 143 galaxies and was analysed 
assuming model DI with Qq = 0.4. Shown is ln(L), and contours are separated by ln(2), 
maximum is indicated by cross and contours are constantly decreasing moving away from 
maximum. 
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Fig. 5. — (0, cz) diagrams of Coma, including data from both zcat and vH. Region surrounded 
by dashed lines is the one selected for maximum-likelihood analysis. See text for details. 
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Fig. 6. — Projected view of Coma. Plotted is (9 cos(0), 0sin(0)) for the same galaxies as 
shown the previous figure. 
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Fig. 7. — Likelihood contours for Coma. Shown is ln(L) and contours are separated by ln(2), 
maximum is indicated by cross and contours are constantly decreasing moving away from 
maximum, (top left) assuming model DI, a cz = 300 km s -1 , and Q = 0.4, (topright) DI 
profile and assuming a cz = 400 km s -1 , and Q = 0.4, (bottum) VI profile and assuming 
a cz = 300 km s -1 . 
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Fig. 8. — Mass estimates at r = 2.5h l Mpc and r turn of Coma based on the DI models. 
Filled symbols at lower radius, circles for a cz = 300 km s -1 , triangles for a cz = 400 km s _1 . 
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Fig. 9. — Mass-to-light ratio at r = 2.5h~ 1 Mpc and r tU m for model DI as function of and 
assuming a cz = 300 km s _1 (filled and open circles) and a cz = 400 km s _1 (filled and open 
triangles). 
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Fig. 10. — Velocity profiles from the models (thin lines) and estimates based on the light 
profile, assuming f2°- 6 /6 - 75 = 1.0, 0.8, 0.6, 0.4, and 0.2 (thick lines, starting from below). 
See text for details. 



